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Spin effects in the phasing of gravitational waves from binaries on eccentric orbits 
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We compute here the spin-orbit and spin-spin couplings needed for an accurate computation of 
the phasing of gravitational waves emitted by comparable-mass binaries on eccentric orbits at the 
second post-Newtonian (PN) order. We use a quasi-Keplerian parametrization of the orbit free 
of divergencies in the zero eccentricity limit. We find that spin-spin couplings induce a residual 
eccentricity for coalescing binaries at 2PN, of the order of 10~ 4 -10~ 3 for supermassive black hole 
binaries in the LISA band. Spin-orbit precession also induces a non-trivial pattern in the evolution 
of the eccentricity, which could help to reduce the errors on the determination of the eccentricity 
and spins in a gravitational wave measurement. 



I. INTRODUCTION 

Up to now most studies concerning gravitational wave 
emission from binaries have been done assuming circular 
orbits (see e.g. jlL an d references therein). Numerical 
studies of the formation of binaries and of their subse- 
quent development suggest instead that the orbits could 
be eccentric 0, |j] ■ It is thus of relevance to take the ec- 
centricity into account when one computes the waveform 
emitted by such a system. This problem has not yet 
been fully explored, mainly due to its complexity, and 
here as a further step we derive the equations which gov- 
ern the time evolution of the orbital parameters, and in 
particular the eccentricity, including the spin-orbit and 
the spin-spin couplings needed for an accurate computa- 
tion in the post-Newtonian approximation up to 2PN of 
the phasing of the gravitational waves emitted by binary 
systems, with components of comparable mass. 

The PN approximation is valid when the two object 
forming the binary are sufficiently separated. The issue of 
estimating the limit of its validity has been tackled with 
different methods: comparisons between post-Newtonian 
templates and results from numerical relativity, for non- 
spinning binaries on quasi-circular orbits |5|; and com- 
parisons between different post-Newtonian waveforms for 
spinning and non-spinning binaries on quasi-circular or- 
bits [g, 0] . All of these studies found a remarkable relia- 
bility of the PN approximation up to separations as small 
as the innermost stable circular orbit, r = QGM/c 2 . It 
is not clear if that still holds for eccentric binaries, and 
this must be answered by extending these comparisons 
to such systems, but one can trust that for low enough 
eccentricities, the post-Newtonian approximation is reli- 
able up to the end of the inspiral. 

We also derive a quasi-Keplerian parametrization of 
the orbit free of divergencies in the zero eccentricity limit, 
and find that spin-spin couplings induce a residual eccen- 
tricity of 2PN order after the orbit has been circularized 
by gravitational wave emission. 

We then solve the equations which govern the evolu- 



tion of the eccentricity and the mean motion for different 
values of the masses and spins as a function of the initial 
eccentricity. 



II. KEPLER EQUATIONS AND EVOLUTION 

OF THE MEAN MOTION AND OF THE 

ECCENTRICITY 

As spin-orbit couplings appear at 1.5PN order and 
spin-spin couplings at 2PN order, it is sufficient to con- 
sider only the Newtonian and spin-coupling term in the 
equations of motion. For simplicity, we will use a sys- 
tem of units where G = c = M = 1, where M is the 
total mass of the system. We start from the generalized 
Lagrangian in the center of mass frame used in [g, |9j : 



£ 



where 
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"Vol ' S2 



(v x a)-£- -5- (as x v) ■ « + £) 



^ (x ■ SO (as • S 2 ) , 



v = m\mi, 
r= \x\, 

C = s l + s 2 , 

mi ?7i2 



The equations of motion are 
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d£_ _ d_ i 

dv l dt 



dp i _ dC 
dt dx l ' 



(1) 



(2) 
(3) 

(4) 

(5) 



(6) 

(7) 



where s l = dC/da 1 . 

We can solve them order by order, which gives at 2PN 
order 
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p = VV + —X X (2£ + £) 



(8) 
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a = — 7T 



x ■ v 



x x (6< + 3£) 



-^vx (4< + 3|) + -J(:r x v) ■ (6C + 6£) 

3t 3 

,S 1 • S 2 , [(x • S 2 ) Si + (a: • Si) S 2 ] 

1 Fit 

+ -^ (x • SO (as ■ S 2 ) . 



(9) 



The reduced energy and reduced orbital angular mo- 
mentum are given by 

J = — (x x p + v x s) 

= x x v H — r-a; x [x x (2£ + £)] — —v x (v x £); 



(10) 



E— -(p-v + s-a — C) 



;S\ ■ S2 



r(x Xtl)'{ 

-^ (x ■ Si) (x • S 2 ) 



(11) 



The magnitude of J is not constant along an orbit [T(| ■ 
Indeed, due to spin-spin interactions, both spin vectors 
undergo a precessional motion and thus, from the con- 
servation of the total angular momentum, it follows that 
J changes at the 2PN order. If we denote its angular av- 
erage (with respect to the true anomaly v, defined later) 
by L, and define A = y/1 + 2EL 2 , we get 



J = L 
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2L 3 
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- 3 cos[2(u -ip)]+A cos(3u - 2^)}, 
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72 



J x Si 



J X S-2 



(12) 
(13) 



where we defined if> to be the angle subtended by the 
bisector of the projections of S, in the plane of motion 
and the periastron line. 

We can find a quasi-Keplerian solution to these equa- 
tions, as (see the appendix) 

r = o(l — e r cosu) + /,-cos[2(w — ip)], (14) 

<j> = (1 + k)v + fa sin(v - 2iP) + / 0)2 sin[2(u - j))], 

(15) 



v = 2 arctan 
I = n(t — t ) 




e t sin u, 



(16) 

(17) 



where (r, 4>) is a polar coordinate system in the plane of 
motion, n is the mean motion, u, v, and I are the eccen- 
tric, true, and mean anomalies, a is the semi-major axis, 



et, e r , and e^ are eccentricities, k accounts for perihelion 
precession, and the f t are constants. 

This paramctrization is different from the one found 
in [llj, which suffered from an apparent singularity in 
the limit e — > (all three eccentricities tend together 
towards zero). This singularity was due to the fact that 
the authors used as a definition of the eccentric anomaly, 
denoted in their paper by £, 

r (0 = 2 [ rmax + r min ~ ( r max - r min ) cos£] , (18) 

which leads to Eq. (ji~4"f with f r = 0. The zero eccentric- 
ity limit of the equations of motion r(<j)) and <f) leads to 
r = f + Srcos[2(cf> — ip)] (see the appendix). If /,. = in 
Eq. (fT4|) . this angular dependence must come from the 
change of variables u(4>). To cancel the e r = 0(e) fac- 
tor in front of cos(u) so that the angular dependence does 
not vanish in the zero eccentricity limit, the function u(4>) 
must be of order 0(e~ 1 ). This is the origin of the ap- 
parent singularity in the quasi-Keplerian parametrization 
found in [ill ]. 

Our paramctrization has the advantage of being free 
from singularities in the zero eccenticity limit, so that the 
latter can be more transparently studied. Note, however, 
that the periastron line (defined by the equation u — v = 
2pir, p € Z) does no longer correspond to r = r mm . 

The mean motion and time eccentricity are given, in 
terms of E and L, as 



n = (-2E) 3 / 2 , 
e 2 = A 2 + —f3 (8, 6 - 2A 2 
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We can invert these relations and find E and L as 
functions of the post- Newtonian parameter x = n 2 ' 3, and 
the eccentricity e = et- These are 
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These allow us to express the constant parameters of 
the quasi-Keplerian motion as 
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We can now use the results from [T(J, [Q , where the 
orbit averages of dE/dt and dL/dt due to the emission 
(26) °f gravitational waves were computed: 
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(40) 
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where ipf is the angle subtended by the projection of Si in the plane of motion and the periastron line. 

We can express these orbit averages in terms of x and e using the post-Newtonian expressions (|23[) and (|24p . Using 



Eqs. (TT9l) and (|20|) . we find the time derivatives of the mean motion and of the eccentricity: 
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We find perfect agreement with [12j, where the spin- 
orbit effects were computed in terms of a and e r . One 
can worry that these derivatives depend on the angles tpi, 
which are not well-defined in the circular limit. This is 
however not a problem, as this dependence disappears in 
this limit both for dn/dt and de 2 jdt. 

We can see that the spin-spin couplings computed here 
induce a positive derivative de 2 jdt for e — > 0. However, 
in symmetrical situations (if the projections of S\/m\ 
and Sijm-2, on the orbital plane coincide), this deriva- 
tive vanishes, due to the fact that r(l, 1,0) — ct(2, 2, 0) = 
(PSi/mi — PS2/m2) 2 , where P is the projection oper- 
ator on the orbital plane. We can compute the value of 
e 2 for which the derivative cancels at 2PN order, which 
is e 2 = 5x 2 [r(l,l,0) - er(2, 2, 0)]/340. We emphasize 
that this effect is independent of the particular quasi- 
Keplerian paramctrization one chooses (see the appendix, 
and in particular Eq. (JA.29I) ). 

We plotted in Fig. Q] the evolution of the eccentricity 
between x — 1/100 and x = 1/6 with spin-orbit and spin- 
spin couplings, for equal-mass binaries with spins uni- 
formly distributed, including also the spin-independent 
PN corrections computed in [l3j , as well as spin-orbit 
precession [14| . We see that spin-orbit precession induce 
a non-trivial pattern in the evolution of the eccentricity, 
which could help to reduce the errors on spin parameters 
in a gravitational wave measurement. We found that the 
quantiles from Fig. [T] are very weakly dependent on the 
mass ratio, whereas the amplitudes of the modulations of 
the evolution of the eccentricity are strongly suppressed 



as the mass ratio decreases. 



CIRCULAR LIMIT 

We define the circular limit of the quasi-Keplerian mo- 
tion discussed above as the limit et —> 0. In this limit, 
we also get e r — > and e^ — > 0. The periastron line 
is not well defined, so that the equations of motion can 
only depend on differences of angles. We find (see the 
appendix) 



x x 



r = x- i + x l ' 2 p{2, 1) + -7i - - 72 cos[2(<?!> - rj>)\ 



%-#•-*«.*,*) 2 
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(44) 

{3 7 i+7 2 cos[2(0-V)]}. 
(45) 



We note that when one includes spin-spin couplings, 
the orbit can no longer be circular in the sense that the 
radius depends explicitly on the angle along the orbit, as 
already mentioned in [8j. This however is not a residual 
eccentricity, as the radius is symmetric with respect to 

4> — > <j) + IT. 

The angular frequency d<j)/dt is not constant. However, 
we can use Eq. (|15[) and define its average along an orbit 
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FIG. 1. Evolution of the eccentricity between x = 1/100 and 
X — 1/6 with spin-orbit and spin-spin couplings, for equal- 
mass binaries, at the top starting from e 2 = 5s; 2 [t(1, 1, 0) — 
<t(2, 2, 0)]/340, and at the bottom from e = 0.01, with spins 
uniformly distributed. In each plot, the grey region is between 
the 5th and the 95th percentile, the solid line is the median, 
and the dashed line is a typical realization. 



as 



t(v=7v) 



t(v — — Tr) 



— dt 
dt 



n 

2^ 

Ti 
— [<j)(v = 7r) -4>{V = -71-)] 



= (-225) 



3/2 



1-^/3(4,3) 



2L 4 



7i 



(46) 



We can thus define a new post-Newtonian parameter 
z = lu 2 ' 3 . In terms of this parameter, the constants E, 
L, and x are 



z z 5 / 2 , , z 3 
£ = - 2 -— 0(4,3)- T 7i, 

_ r 3/2 
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(47) 
(48) 

(49) 



Now, we can use Eqs. 06J, (|32J), and (|33J to find 
dw 961/z 11 / 2 



dt 



3/2 

1-^-/3(113,75) 



2 2 

^ f r(247,721,0) + ^r(7,l,0) 



(50) 



which is in agreement with what was previously com- 
puted in [9|, ll5j . 

Alternatively, if we define a circular orbit to have 
de 2 /dt = 0, which implies e 2 = 5z 2 [r(l, 1, 0) - 
ct(2,2,0)]/340, we get 



du _ 96UZ 11 / 2 

~dt ~ 5 
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1216 



cr(6519, 18 527,0) 



2432 



r(439, 287,0) 
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III. CONCLUSION 

The main result of this paper is the derivation of the 
spin-spin effects in the evolution of the mean motion and 
of the eccentricity, for binaries with an arbitrary eccentic- 
ity. Particularly, the fact that spin-spin couplings may 
induce a residual eccentricity can be important for pa- 
rameter estimation when gravitational wave observations 
will be made possible. If eccentric templates allow to 
measure eccentricities of O(10 _3 -10~ 4 ), the modulation 
induced by spin-orbit precession could significantly im- 
prove the determination of the spins of the binary. 

We also derived the equations of motion r(t) and 4>(i), 
for black hole binaries of comparable mass with Newto- 
nian, spin-orbit, and spin-spin terms on eccentric orbits, 
and found a family of parametrizations free of divergen- 
cies in the circular limit e — ► 0. 
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Appendix: A two-parameter family of 
quasi-Keplerian parametrizations 



We can find a two parameter family of quasi-Keplerian 
parametrization of the equations of motion by taking the 



ansatz 



r = a (1 — e r cos u) + f r i cos(i; — 2ip) 

+ / r , 2 cos(2y-2V>), (A.l) 

<j>=(l + k)v + U,i sin(v - 2ip) + U i2 sin(2w - 2ip) 

+ U, 3 sm(3v-2iP), (A.2) 



v = 2 arctan 



±±**tanH 

1 - e A 2 



(A.3) 



I = n{t — to) = u — e t smu + ft sin(v — 2ip), (A. 4) 

where (r, 0) is a polar coordinate system in the plane of 
motion, n is the mean motion, u, v, and I are the eccen- 
tric, true, and mean anomalies, a is the semi-major axis, 
et, e r , and e^ are eccentricities, k accounts for perihelion 
precession, and the fi are constants. 
The constants are 
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(A.10) 
(A.ll) 
(A.12) 
(A.13) 
(A.14) 
(A.15) 
(A.16) 



where Ai and A2 are arbitrary functions of A. 

The Kcresztes-Mikoczi-Gergely (KMG) quasi- 
Keplerian parametrization of the orbit [llj is obtained 
by imposing f r j = 0, which leads to A2 = and 
Ai = -1/A 

To get Eqs. (fT4"f to (IT7|) . we required that the values 
of the semi-major axis and of the eccentricities should 



not depend on the position of the periastron line in the 
orbital plane. This implies Ai = A2 = 0. 

We can express E and L as functions of the post- 
Newtonian parameter x = n 2 ' 3 and of the eccentricity 
e = e t . We get 
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The constants in the equations of motion then become 
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One can see that for the solution to be free of diver- 
gencies in the zero eccentricity limit, Ai and A2 must be 
regular as e — > 0, which implies, at this PN order, that 
they must be regular functions of A. It is not the case for 
the KMG parametrization, but this is in fact a coordi- 
nate singularity, since, as we will see later, the equations 
of motion have a well-defined zero-eccentricity limit. 

The effect of Ai and A2 on the evolution of the eccen- 



tricity with respect to Eq. (J43I) is 



leading order in e, 



.de 2 



' y - r cr(0,0,688e + 2139e 3 + 148e 5 ), 



* 15(l-e 2 ) 9/2 

(A.29) 
which, as long as the parametrization is regular, does not 
affect the residual eccentricity we found in this paper. 



CIRCULAR LIMIT 



Let us now compute the equations of motion r(4>) and 
d<j>/dt. From Eq. (|A.2|) . we get 



v = (1 - k)(j) - fax sin(</> - 2i/>) - f<p,2 sin(2</> - 2ip) 
-U, 3 sm(3<j>-2ip). (A.30) 



Then, together with Eq. (|A.3j) . Eq. (|A.lj) becomes, at 



r = x- 1 + x^ 2 /3(2, 1) + | 7 i - |72{ cos[2(0 - 4,)] 

+ (A 2 + e 2 Ai) cos(</> - 2^) - e 2 A x cos(0 + 2^)}, 

(A.31) 

and we can compute from Eq. (JA.4I) . at leading order in 



d<f> 



3x 7 / 2 



x 3 / 2 -x 3 /3(4,3)-^-7i 
+ 12X 7/2 [ (A 2 + 2e 2 Ax) cos(0 - 
- 2e 2 Ai cos(> + 2ip)] . 



.7/2 



2 
2V) 



-72 cos[2(0 - ip)] 



(A.32) 



One can see that the use of the KMG parametriza- 
tion does not change the equations of motion in the zero 
eccentricity limit with respect to the one obtained from 
Eqs. (fT4|) to (fl7|) . Furthermore, as the periastron line is 
not well-defined in the circular limit, we have to impose 
that A2 — > as e — > (which is equivalent to A2 ~ 0(A)), 
so that the equations of motion in the circular limit are 
independent of the choice of an arbitrary periastron line. 
From Eq. (IA.2[) . we see that such an arbitrary choice in- 
duces a non-zero value for 4>(v — 0) of 2PN order. How- 
ever, the differences induced in the equations of motion 
are subsequently of 4PN order, far beyond the limit of 
their validity. 
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